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Dynamic nuclear polarization (DNP) is a promising strategy for generating a significantly in¬ 
creased non-thermal spin polarization in nuclear magnetic resonance (NMR) applications thereby 
circumventing the need for strong magnetic helds. Although much explored in recent experiments, 
a detailed theoretical understanding of the precise mechanism behind DNP is so far lacking. We 
address this issue by theoretically investigating solid effect DNP in a system where a single electron 
is coupled to an ensemble of interacting nuclei and which can be microscopically modelled by a 
quantum master equation. By deriving effective equations of motion that govern the polarization 
dynamics we show analytically that DNP can be understood as kinetically constrained spin diffu¬ 
sion. On the one hand this approach provides analytical insights into the mechanism and timescales 
underlying DNP. On the other hand it permits the numerical study of large ensembles which are 
typically intractable from the perspective of a quantum master equation. This paves the way for a 
detailed exploration of DNP dynamics which might form the basis for future NMR applications. 


PACS numbers: Valid PACS appear here 


The sensitivity of imaging and spectroscopy techniques 
based on nuclear magnetic resonance (NMR) depends on 
the spin polarization that arises from the Zeeman inter¬ 
action of the nuclei with an externally applied magnetic 
field. However, even with the highest magnetic fields 
that can currently be generated by superconductive mag¬ 
nets, the nuclear spin polarization of nuclei is only in 
the order of less than ten parts per million at ambient 
temperature. A promising strategy to overcome the low 
sensitivity in many NMR applications is the use of dy¬ 
namic nuclear polarization (DNP). Here the electronic 
spin polarization — that is about three orders of magni¬ 
tude higher than the nuclear spin polarization — is trans¬ 
ferred from paramagnetic centres in solid state materials 
to the surrounding nuclear spin ensemble by irradiating 
with microwave fields at an appropriate frequency. DNP 
was already proposed and demonstrated in the early days 
of nuclear magnetic resonance [11-0, but it has recently 
attracted increased attention since progress in hardware 
technology, mainly in the form of stable high frequency 
and high power microwave sources, has made it possible 
to substantially enhance the NMR signal in several appli¬ 
cations Hi including the detection of nuclei with small 
gyromagnetic constants 0, of nuclei located on the sur¬ 
face of structured materials 0 and of nuclei in protein 
micro-crystallites 0 by magic angle spinning solid state 
NMR spectroscopy. Furthermore, DNP carried out at 
cryogenic temperatures combined with a fast dissolution 
step to rapidly bring the sample to ambient tempera¬ 
ture, was successfully used to generate highly polarised 
^^C-labelled compounds that could be detected in in vivo 
Magnetic Resonance Imaging applications after their in¬ 
jection ids. A full understanding of the spin dynam¬ 
ics during DNP is essential for the interpretation of the 
experimental data and ultimately for optimising the ex¬ 
perimental protocols. In principle this problem can be 


addressed by numerically solving the full quantum mas¬ 
ter equation Owing to its complexit y, h owever, 

this approach is limited to small systems 12|, [l31 . 


In this work we show how to overcome this limita¬ 
tion and to extend the analysis of the spin dynamics to 
large ensembles containing more than a thousand cou¬ 
pled spins. By using an adiabatic elimination procedure 
we derive an effective equation for the evolution of the 
spin polarization. This equation is classical in the sense 
that it describes transitions between Zeeman basis states 
and therefore the computational complexity of simulat¬ 
ing the spin dynamics is drastically reduced. Moreover, 
the analytical form of the effective equations of motion 
permits a transparent interpretation of the physics un¬ 
derlying DNP in terms of operator-valued, i.e kinetically 
constrained, diffusion rates. 


Solid effect DNP (SE DNP) is observed in systems that 
consist of interacting unpaired electrons S and nuclei / 
[see Fig. [TJa)] in an external magnetic field with Lamor 
frequencies ws and uJi. Polarization buildup is mediated 
by the electrons that are irradiated at an off-resonance 
frequency coq which is in resonance with either the so- 
called double or zero quantum transition ws ± 14 - ^ . 
For SE DNP to be effective, the electrons need to have a 
resonance linewidth that is smaller than the nuclear Lar- 
mor frequency. A suitable minimal model for SE DNP 
consists of an ensemble of n nuclei coupled to a single 
unpaired electron through the hyperfine interaction. The 
spin dynamics of SE DNP can be described by the quan¬ 
tum master equation dp/dt = Lp with the Liouvillian 


L = —iH -l-F. It describes the coherent evolution of the 
spin system under the action of the Hamiltonian commu¬ 


tation superoperator H and the effect of the dissipative 
processes (decoherence) represented by the relaxation su¬ 
peroperator r. Note, that throughout a bold symbol is 
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FIG. 1. (a) Sketch of the solid effect DNP setup. The elec¬ 
tron is coherently driven by the applied microwave field and 
coherent excitation-exchange between the spins is mediated 
through dipole-dipole couplings. Furthermore the electron 
and nuclei are subject to dissipation leading to relaxation 
(see text for details), (b) The effective dynamics of the polar¬ 
ization is described by a classical master equation and thus 
permits the simulations of large systems, (c) Comparison be¬ 
tween the exact (solid black lines) and effective (colored lines) 
relaxation dynamics of a system consisting of one electron and 
four nuclei, (d) Polarization dynamics in a linear spin chain 
with one electron and 30 nuclear spins. The polarization ver¬ 
sus spin index (0 denotes the electron position) is shown for 
different times and the polarization versus time is shown for 
the electron (blue line) and four selected nuclear spins (red 
lines). 


used to represent superoperators and a doublehat for a 
commutation superoperator, i.e. O = [O, •]. 

In the following we discuss the structure of the quan¬ 
tum master equation in detail. The Hamiltonian H of the 
system consists of four terms: H = + Ho + H+ -|- 

where Hi = uJi{Sz J2k^kz) is the Zeeman term. The 
term 

Hq = xSz + AkikzSz + dkj (^ikzijz — Ik ■ ij^ 

k k<.j 

describes the hyperfine and nuclear dipole interaction. It 
depends on the secular strengths of the electron-nuclear 
hyperfine interaction, parameterized by Ak and further¬ 
more on the inter-nuclear dipolar coupling coefficients djk 
[^ . There is an additional dependence on the so-called 
offset parameter A = ws — wi — wq which represents a res¬ 
onance condition between the microwave frequency wq 
and that of the double quantum transition ujs — uji. The 
remaining terms of the system Hamiltonian are 

k 


which arise from the pseudosecular hyperfine interaction 
and the microwave irradiation. They are parameterized 
by the pseudosecular coupling strength Bk± between the 
electron and the nuclei and the microwave amplitude Wi. 
Note, that the interaction parameters Ak, Bk+ and dkj 
depend on the geometry of the spin system [21| . 

The relaxation superoperator is composed of two 
terms, T = Ti -|- r 2 . The first one is given by 

ri = ^[D(5+) + D(5_)] + rth 

d(I) 


where we have defined the dissipator D(X)/3 = XpX^ — 
{p,X^X}/2. Fi describes dissipative processes or lon¬ 
gitudinal relaxation due to interaction of the spins with 
the environment 22|. The rates for longitudinal relax¬ 


ation are and i?® for the electron and the nuclei, 

/g\ 

respectively. The superoperator Fth = Hoi?) '[D(5'_) — 
T>{S+)]/2 acts as a correction term and ensures that 
the electron relaxes back to its thermal equilibrium state 
which—at typical temperatures of T > O.IK— is up to 
a scaling factor approximated by pth = 1 — ‘IPqSz- Here 
Pq = tanh (/kas/2fcBT) is the thermal electronic polar¬ 
ization [^. The second term. 


r2 = 2i?f D(5,) + 2E4'^D(4,) 

k 


represents the loss of coherences or transverse relaxation 
of the spin ensemble, which depends on the rates i?) and 
i?® for the electron and the nuclear spins, respectively. 

The full quantum master equation describes the 
dynamics of the spin system in the Liouville space 
L of dimension 2^^ where N = n -I- 1. How¬ 
ever, all information that is relevant for character¬ 
izing SE DNP is typically contained in the Zeeman 
subspace £z [1^ which is spanned by the operators 
{i, ikz, ikzik'z, Ikzlk'zh"z, ■ ■ ■} and has a dimension 
2^. For example the bulk nuclear polarization is given 
by the expectation value of Ikz- We have obtained a 
closed equation for the dynamics in the Zeeman subspace 
by means of a quite lengthy adiabatic approximation pro¬ 
cedure which is detailed in the supplemental material. 
This resulting effective master equation is of the form 
dp 2 ,/dt = Lz/5z. It has not only a substantially reduced 
dimension but in addition it is of entirely classical nature 
— although we employ in the following a notation that is 
also used in the quantum case. As shown in the supple¬ 
mental material the dynamics can be expressed in terms 
of jump operators [see Fig. nib)], which connect classi¬ 
cal spin configurations, and their corresponding effective 
jump rates, i.e. the equation is of Lindblad form [ 2 ^. 
Similar approaches have been employed in the study of 
cold atomic systems When formulated in this 
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way the problem lends itself to the use of classical ki¬ 
netic Monte Carlo (kMC) algorithms which permits the 
simulation of large systems. To check the validity of this 
approximate approach we have simulated small systems 
with the exact quantum Master equation and kMC. Ex¬ 
ample results are shown in Fig. [Ijc) as well as in the 
supplemental material and show excellent agreement be¬ 
tween the two methods. 

The effective dynamics in the Zeeman subspace is gov¬ 
erned by a set of single-spin and two-spin dissipators: 
Lz = Lsingie-spin + Ltwo-spin- The single-spin terms are 
given by 


-^single—spin 


= rv^D(5+) -f rL^D(5_) -f 


-Et 


r«D(4+) + r;,^^D(4_: 


( 1 ) 


and depend on the constant rates 


r^s) 


p(I) _ 

^ k± — 


1 T Po p(S) 

r, Pi 


H 


(I) 


Bh B, 


r(S) 

k+ p(I) 

- JXry 




( 2 ) 

(3) 


fs) 

Here T^ are the rates at which the electron spin flips 

’up’ (-I-) or ’down’ (—). They depend on the longitudi- 
(S) 

nal relaxation rate R\ of the electron with a weighting 
pre-factor which ensures that in the absence of any per¬ 
turbation by the microwave field the electronic steady- 
states polarization is identical to the thermal electronic 
polarization Pq. The second term is a consequence of the 
applied microwave field. Eq. © describe the rates for 
a nuclear spin k to flip ’up’ or ’down’. Since we have 
assumed in the derivation of these equations that the 
thermal nuclear polarization is negligibly small, these two 
rates coincide and correspond to half the longitudinal nu¬ 
clear relaxation rate i?®. The additional term appearing 
in Eq. ([3|) arises from the pseudosecular interaction with 
the electron which leads to a tilt of the nuclear effective 
field axis. 

The two-spin dissipators are given as 


L+wn — 


two—spin 


= (4) 


fc<j 


with Yfc = h+S- -I- Ik-S+ and Xkj = h+Ij- + Ik-Ij+ ■ 
They clearly describe a diffusive (flip-flop) dynamics be¬ 
tween the electron and the nuclei as well as among the 
nuclei. Interestingly the diffusion rates T^^®^ and f 
are not constant, but operator-valued. Hence the diffu¬ 
sion rate of the spin-polarization depends on the state of 
the entire spin system. In the context of the study of 
glassy relaxation 2^, ^ such operator-valued rates are 
often referred to as kinetic constraint. Explicitly, they 


read 


p(IS) _ 

^ k — 


Dk = 


\Bk-B, 


k+\ 




Al -I- ^slsz 




-1 


(5) 


p(ii) _ 


m. 


(I) 


(i + dt) 


Ckj — 


{Ak - A,)S, 


2R. 


(I) 


( 6 ) 


A brief discussion of these expressions provides insight 
into the spin dynamics on which SE DNP is based on. Eq. 
dD) describes the rate of flip-flop transitions or jumps be¬ 
tween a nuclear spin k and the unpaired electron. Hence 
r), is of fundamental importance for the SE DNP effect. 
This rate depends on both the strength of the microwave 
field (jJi and the pseudosecular interaction strength Bk± ■ 
It is inversely proportional to the square of the nuclear 
Larmor frequency wi and also inversely proportional to 
the sum of the transverse relaxation rates of the electron, 
R^\ and the nuclear spin, R^^. This indicates that al¬ 
though this process is represented by the flip-flop jump 
operator, the underlying quantum mechanical process is 
mediated by coherences which in fact decay with these 
relaxation rates. 

Most importantly, the rate for polarization transfer be¬ 
tween the electron and the fc-th nucleus is controlled by 
the operator Dk whose (eigen)value depend on the polar¬ 
ization level of the nuclear spin ensemble. For negligible 
nuclear polarization in an ensemble of nuclear spins the 
sum J2s^k ^shz evaluates to a value close to zero as 
there is an equal probability for the nuclei to be in their 
spin ’up’ or ’down’ state. However, with increasing level 
of nuclear polarization the magnitude of this sum grows 
which decreases the effective rate T), . Note, that the ef¬ 

fect of the off-resonance parameter A is negligible as long 
as the microwave frequency wq is matched to the double 
quantum transition at ws — wj. 

The transport of polarization within the dipolar- 
coupled network of the nuclear spins takes places via 
spin diffusion which is governed by the operator-valued 
rate given in Eq. ([HI). The rate of polarization trans¬ 
port among two nuclei is proportional to the square of 
their nuclear dipolar coupling constant dkj, and it is in¬ 
directly proportional to the nuclear transverse relaxation 
rate i?®. Interestingly, there is also a dependence on the 
difference of the secular hyperfine interaction strength 
Ak of the two nuclei. In the case in which this difference 
is significant, the internuclear polarization transfer be¬ 
tween the spins can become strongly quenched provided 
the electronic polarization Sz is non-zero. 

As the effective dynamics in the Zeeman subspace can 
be treated efficiently numerically SE DNP can be stud¬ 
ied in relatively large samples. In Fig. [TJd) we already 
showed a simple example for 31 spins in one dimension 
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FIG. 2. Simulation of the polarization dynamics of 1330 nuclear spins (^^C) arranged on a regular grid in a cube (11 x 11 x 11) 
around one central electron (black). The diameter and the colour of the spheres indicate the expectation value of the nuclear 
spin polarization at the given position. The simulation demonstrates that polarization buildup is fastest for the nuclei that are 
located close to the electron at an angle of 7r/4 with respect to direction of the applied static magnetic field. In order to show 
the contribution of spin diffusion to the distribution of polarization within the spin ensemble the dipolar coupling constants 
were set to zero in the rightmost panel. This prevents any transport of polarization through spin flip-flop transitions between 
pairs of nuclear spins. Note, that the front quarter of the spin ensemble is not shown for better visibility of the polarization 
buildup near the electron. 


while Fig. [2] now shows the progressive buildup of spin 
polarization within a three-dimensional sample of 1330 
spins arranged on a regular cubic lattice with one 
electron spin at the centre. For the latter the distance 
between the nuclei is 10 ± 5% A and the field strength of 
the static external magnetic Bq field is 3.4 T. The applied 
microwave field has a frequency of wq = ws — = 34.5 

MHz and a strength uji =100 kHz. The relaxation 
time constants are = 1 = 10® = 

1.4 X 10-4 s-4 and = 10^ s”!. 

Let us finally investigate the role the nuclear dipolar in¬ 
teraction plays in the distribution of polarization between 
the nuclear spins. This ques tion has been addressed e.g. 
by Hovav et al. in Ref. [ll| who conducted simulations 
of small chains (up to 9 nuclei) within the framework of 
the full quantum master equation. Their results show 
that there is only a weak dependence of the polarization 
build up rates on the nuclear dipolar couplings between 
nuclei further away from the electron (nuclei belonging 
to the bulk) suggesting that there is a direct transfer of 
polarization from the electron to all bulk nuclei through 
mixing of the multi-nuclear product states. Our data 
shown in Fig. [5] clearly suggests that the nuclear dipo¬ 
lar couplings are important. To further investigate their 
role in polarization transport we have simulated the po¬ 
larization dynamics of a linear spin chain consisting of 
30 spins as shown in Fig. [TJd) [see supplemental ma¬ 
terial for further details]. In contrast to the analysis in 
Ref. [m our findings show that a decrease of the nuclear 
dipolar interaction constants of the bulk nuclei by a fac¬ 
tor of two alters the polarization buildup dramatically. 
This is consistent with a diffusive polarization transport. 
The discrepancy is thus most likely caused by boundary 
effects, which can be severe for short chains. Moreover, 
further discrepancies might arise due to the way in which 
relaxation was introduced in the reference frame in Ref. 


Ill] as was already previously discussed in Ref. [22l| . 

Let us now establish a more quantitative connection 
between SE DNP and the diffusion of polarization. For 
a linear chain of n nuclear spins, as shown in Fig. [l] we 
can derive an average diffusion constant: 


n — 1 




= - 7 X! ^ 

r? — 1 


^R2^dlal 




^4i?2 + {Ak — Ak+i) 


( 7 ) 

Here is the distance between the nuclear spins k and 
fc -|- 1, dk is the dipolar coupling constant between the 
two spins and Ak is the strength of the secular hyperfine 
interaction of spin k. The constant Dk is the rate of dif¬ 
fusive polarization transport among neighboring nuclei. 
It is quenched when the difference between the secular 
hyperfine interactions of the two spins, \Ak — ^fe-i-i|, is 
large and when the nuclear transverse relaxation rate R^'^ 
is fast. Comparing the polarization dynamics calculated 

q2 

via kMC with a simple diffusion model, dp/dt = 
where p = p{x, t) is the polarization field, indeed shows 
good agreement with data as the one presented in Fig. 
[Ud). Further details and a discussion of the relevant 
boundary conditions is provided in the supplemental ma¬ 
terial. We conclude that polarization transport in a SE 
DNP experiment is mediated by nuclear spin diffusion. 
The presence of the electron will hinder — or constrain 
— this process for nuclei close to the electron since for 
these nuclei the difference between the strengths of the 
secular hyperfine interactions is significant IT, 3^. 

We have derived effective equations of motions which 
show that the mechanism underlying SE DNP is kineti- 
cally constrained diffusion. The derived set of equations 
enable the study polarization buildup in relatively large 
spin systems which has been demonstrated by investigat¬ 
ing DNP in a system with 1330 nuclei and a single 
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electron. Finally, we have highlighted the importance of 
nuclear dipolar couplings in SE DNP and derived an ef¬ 
fective rate for polarization diffusion in a one-dimensional 
chain. We believe that this study represents an impor¬ 
tant step in the quantitative understanding of SE DNP 
for practical applications in NMR. In the future we will 
investigate whether the method presented here is also 
applicable in more complex settings, which involve e.g. 
more than a single electron. 
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SUPPLEMENTARY MATERIAL 
Adiabatic Elimination 

For any linear system with constant coefficients 

(7 = La (8) 

defined in some vector space "H, the following technique is valid. 

Suppose we need to know dynamics in a subspace "Hi C Ti and initially also (t(0) = (Ti(0) G Hi- By using the space 
decomposition 


n = ni + n2 


and introducing the projections 


dfc — ^kj — ^kljTTj, kjj — 1, 2, 


we can write 


a\ — LiiCTi + Li2d'2, d'2 — L22d'2 + L2lfTi. 


Given a \, the second equation is resolved as 

(T 2 (t) = f e“^^^'^L 2 i(Ti(T) dr. 

Jo 

Substituting this into the first equation, we obtain the equation closed in "Hi 

(Ti(t) = Lii(Ti(t) + / K(t - r)(Ti(T)dr, 


where 


(9) 


K(T) = Li2e^-^L2i. 

In terms of the Laplace transforms (defined for C > 0) 


^+oo I'+ca 

hiC) = / e~‘^*ai(t)dt, 1k(C) = / e~^*K(t) dt, 

Jo Jo 


is resolved as 


where 


/i(C) = (P(C))-'di(0), 

P(C) = C1-Lii-1k(C). 


( 10 ) 


( 11 ) 


( 12 ) 


The initial value (ri(0) and the Laplace transform of the kernel Ik uniquely define the Laplace transform of the 
solution. The solution itself is then uniquely found by inversion of its Laplace transform. 

It follows from ® that the solution at time t depends on its history in the time interval [0,t]. The integral kernel 
K is often called memory function. Besides the memory effects, K describes how dynamics in the complementary 
subspace 'H 2 affect the dynamics in Hi- 

If all eigenvalues (k of the superoperator L 22 have negative real parts then any solution to ® tends to the steady- 
state solution 


di{t) —y t —>■ + 00 , 


satisfying 


(Lii -|- lV[)(Ti^. — 0 


7 


with 


p+OO 

M= K(t)dt = -Li2L22”^L21. 

Jo 

On the other hand, in terms of (fT0|l . (fT2|l and the Laplace transform, 

K(t)=^eC'=‘Kk, 1 k(C) = ^^Ll_. 

Due to the formulas 

W0) = -Ef 

we have that under the conditions 

ICI<C-, IIL12II ■ IIL21II < C-, C-=min|Cfe| 

we obtain the approximation 

1k(C) = 1k(0) + ^(0)C + ...^M 

where the second term is negligibly small. Since the small values of ^ are responsible for the long-term slow dynamics, 
it means that the steady-state part M of the integral term well approximates the informative long-term asymptotic 
of the dynamics. 

Thus, under the conditions 


||Li2M|L2i||«C-, ||Lii+M||«C-, 


(13) 


we have the asymptotics in (1121) 


Returning to the Laplace transform dl), this means that solutions to equation ([9]) are well approximated by solutions 
to the equation 


(Ti(t) = (Lii -f M)(Ti(t). 


(14) 


Equation (fT4l) does not contain any memory function. Its operator Ln -|- M is time-independent. Due to (fT^ . fast 
dynamics in the complementary subspace 7^2 make the solution in Tii tend to the steady-state, rapidly “forgetting” 
its previous history. 

We call equation (1141) an adiabatic approximation for equation ([5]) and say that the complementary subspace 772 
is adiabatically eliminated. The term adiabatic means that we first separate fast motions in 'H 2 from slow motions 
in 77i and then eliminate all information about the fast motions, not influencing the slow dynamics. Here the 
equation preserves its initial form of the linear system with constant coefficients. The technique is easily generalised 
to inhomogeneous systems and initial conditions outside the informative subspace 77 1 . 

Such procedures are commonly used (among other fields) in condensed matter theory, quantum optics (where the 
term adiabatic elimination initially comes from), quantum statistics and quantum information, in cases where either 
external driving or an interaction with an environment are involved in such way that a non-resonant part of the 
dynamics becomes non-informative. A similar effect, well-known in the NMR context, is achieved by proceeding to 
a rotating frame and averaging out secular terms. Another example is eliminating high spin correlation orders in 
multi-spin dynamics where the role of the fast motion is played by spin relaxation. See in this context, for example, 
ilElIH. If the dimension of the informative subspace 77i is much smaller than the dimension of the total space 
77, this gives a significant reduction of the initial problem. 





Elimination of non-zero quantum coherences 


The Liouville state space can be decomposed as 

where Hq is the subspace of g-quantum coherences, 

H zp = q<^ip, p€Cq, (? = 0, ±1, ..., ±(n + 1). (15) 

Within the method described in the previous section, let the subspace Hi be the subspace of zero-quantum coherences 
Co and H 2 the complementary subspace, 

ni=Co, H2=Y.Cq. 

q^O 

Starting from the thermal equilibrium, we have 

di(0) = 1-2 Po5,. 

The relations 

{iH 0 -I- r)£q C Cq, H ±Cq C Cq±l 

along with (ITSl) give 

Liidi = -{iH 0 -I- r)(Ti, L 211 T 1 = -i{H + + H -)(Ti, 

Li 2 '^±i = —iH =p£±i, \-112Cq = 0, |g| > 1, , , 

^■2iC±i = —(iiw/l -|- iH 0 -I- iH ± + r)£±i, 

^iiCq = —{iquji\ -\- iH 0 -|- iH _|_ -|- iH _ -|- r)/3g, jgl > 1. 

Due to the presence of relaxation, the eigenvalues of the superoperator L 22 have negative real parts, satisfying the 
requirement outlined in the previous section. Physically, at high magnetic field, 

llwf^rll ^e«l, (17) 

so the superoperator M is found as the series 

IVE = UJj ^IVIi -t- U)j ^]Vl2 -b ... (fS) 

converging exponentially fast 

m 

k^l 

Using CS), we can calculate 

Mi = -i[H +,H _], M2 = -H +{iHo + T)H _- H _{iHo + T)H +, ... (19) 

On the other hand, it follows from (TlBl) 

C- > |w/|, IILi2 II < max II 17 ±11, ||L 2 i || < 2 max ||17 ±||. 

Hence, because of dUl), (HID, (HU), conditions m are satisfied, so we can replace equation m by its adiabatic 
approximation 

d = Lod (20) 

with 

Lq = Lii -b IVI = IVlo “b toj h/Ei “b toj- 1V[2, 

where 

Mq = —{iH 0 -b r). 

Ml .2 are defined in (HU), and in (1181) we restricted to the second order approximation for M. 
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Elimination of non-Zeeman spin orders 

In the previous section, we adiabatically proceeded from the initial Liouvillian L to the new Liouvillian Lq, well 
describing the dynamics closed in the subspace £o of zero-quantum coherences. 

We decompose now the zero-quantum subspace as 

-Co = C-z + -Cc, 


where 


— SpClTl^l^ Ikz-) ^kzlk'z-, ^kz^k'z^k"zi • ■ •} 

is the subspace of Zeeman spin orders, and Cc is the complementary subspace consisting of non-Zeeman zero-quantum 
coherences. 

The commutation character of the notation O = [O, ■] implies that the superoperator Mi is a commutation. 

Hip = + g Bk+Bk-Ikz, H'l = —— '^^{Bk+Ik+S- + Bk-ik-S+). 

We have also 


Ho = Hop + H'o, r = ro + r'. 


Hop — XSz + AkikzSz, 


H'o=Y.<^. 

j<k 


kj 


li 








kz 1 


p/ _ ^1 


(5 +1 _ +1 +) + i ^ r[^'> (I kJ k-+I k-I fc+) + Tth. 


( 21 ) 


The superoperators H o_o, H 1 , 0 , Tq trivially act on Cz- Within the accuracy of ~ lOj only the superoperators 
H Q, H I contribute to L 21 , L 12 . The superoperator 

X = —iH 0,0 ~ To 

maps Cc to itself with eigenvalues C'k satisfying 

Re a < 0, IC^I > min{2Ri^\ ^ 

Thus, under the conditions 

mi„ { (24'>)^ (Br> + {t|!, , («<'>)“}, ( 22 ) 

conditions (IT^ are satisfied and the subspace Cc is adiabatically eliminated. We come then to the adiabatic approx¬ 
imation 


(T — Lzct, Lz — Lo ll + — ~Lo,i2Lo,22 ^Lq^ 


(23) 
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closed in the Zeeman subspace Hz- In accordance with the previous notations, 

Lo.kj = k, j = 1, 2, 

where tt^ 2 are projections onto the subspaces Cz,c respectively. 

Using conditions (Ull) , ([22]), it can be shown that the right-hand side of (l23ll is well approximated as 

Lo,ii = -r'-r", 

j<k k 


where F' is given by (EU, 

1 


F" = 




d(U 

) (5 +5 _ + 5 +) + ^ ^ Bk+Bk- (I kj k-+I k-I k+) 


and 


x'„^ = ^{h+ij- + ik-ij+), 


V' — 

I h — 


c'., = 


2R 


0^1 

4w/ 

(U 


iBk+ik+S-+Bk-ik-S+), 


R 




(24^^+ Cl/ 


Ckj = {Ak - A,)S, + —- B,+Bj_)l 

OU)j 

Dk = Ai -h Agisz + (4^1 ~ Bk+Bk-) i. 


s^k 


(25) 


The advantage of formulas is that they reduce the inversion ho 22 ^ in the subspace Cq to the much simpler 
problem of inversions C'f.^, D'j, of Zeeman operators. The latter are defined in the 2^+^-dimensional Hilbert space and 
are diagonal in the usual Zeeman basis. 

The last terms in the expressions for Ckj, Dk are the second order corrections with respect to the inverse nuclear 
Larmor frequency . For typical random spin geometries, these terms are either quenched by the first order terms 
or small simultaneously with them. Hence, these terms can often be neglected. 


The Lindblad form 


Introducing the notation 


and utilizing the double-commutator character of approximation (1241) as well as the thermal correction Fth, it is 
straightforward to see that the right-hand side of the Zeeman projected equation (1231) is written in the purely Lindblad 
form 


Lz = Ff D(5+) + (rSD(4+) + hi^l-Dilk-)) + 

k 


k k<.j 
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with the constant rates 


T.(S) iT-Pon(S) , ^(s) 

U -^—^1 > 


— If j_^ 

ife ±- 2 1^1 + 4^2 -2 


fcl ^(/) 


related to the single-spin jump operators S±, Ik±, and the operator rates 

■^iBk+Bk- 


f (/S) _ 
^ fe — 






pi-* ^ ) _ 


dl, 


ARX 




-Dfe = 


+ Ss/fc "4s/, 




Ckj — 


{Ak - Aj)S, 


2R. 


(I) 


related to the two-spin jump operators 

hfc = Ik+S— + ik-S+, Xkj = ik+ij- -\- ik-Ij+- 


Comparison between full Liouville space simulations and kMC simulations 

We have carried out a number of comparisons between simulations of the polarization dynamics using the quantum 
master equation and our proposed method that involves an adiabatic approximation procedure to obtain a master 
equation for the states contained in the Zeeman subspace Cz and the use of a kinetic Monte Carlo (kMC) algorithm 
to calculate the polarization dynamics. As long as the parameters are chosen in accordance with the conditions 
described in the supplementary material, an excellent agreement between the two simulation methods is achieved. An 
example of such a comparison is given in Fig. ([3]) which shows the polarization dynamics in a system containing one 
electron and four nuclear spins forming a coupled network. The agreement can be further improved by increasing the 
number of averaged trajectories for the kMC method, however this is done on the expense of further increasing the 
computational time. The parameters used for the simulation are summarised in the following table. 


Parameter 

Value 

Nuclei 


Bo 

3.4 T 

Temperature 

1 K 

Number of trajectories 

W 

CJI 

50 kHZ 

nuclear X-coord. 

3.51, -2.52, 5.98 -4.01 (A) 

nuclear Y-coord. 

0.02, 0.03, 0.00 0.00 (A) 

nuclear Z-coord. 

3.56, -4.34, -0.52 0.02 (A) 

electron coord. 

-0.03, -0.01, -0.01 (A) 

\B+\ 

0.935 MHz, 0.822 MHz, 92.1 kHz, 22.2 kHz 

A 

0.318 MHz, 0.794 MHz, -0.352 MHz, -1.26 MHz 

di,2, di,3, di,4 

-0.055 kHz, -0.660 kHz, 0.048 kHz 

^2,3, ^2,4, <^3,4 

0.037 kHz, -1.040 kHz, 0.059 kHz 

Tie, T2e, tin, t2n 

1 s, 10 /rs, 1 h, 5 ms 
























12 



FIG. 3. Comparison between a simulation of the polarization dynamics during SE DNP using the quantum master equation 
in the full Liouville space and a simulation based on a master equation for states contained in the Zeeman subspace that was 
obtained by adiabatic elimination and solved by kMC. For the kMC simulation 10^ trajectories were averaged. The polarization 
dynamics is calculated for a system containing one electron (blue) and four coupled nulcei. There is excellent agreement between 
both simulations with the coloured lines indicating the results obtained using kMC and the black lines showing the results using 
the quantum master equation. The agreement can be further improved by increasing the number of averaged trajectories. 


SE DNP and Spin Diffusion 


To evaluate the role that the nuclear dipolar interaction plays in the distribution of the polarization within a 
larger ensemble of coupled nuclear spins, we have simulated the polarization buildup dynamics in a linear chain of 40 
coupled spins with an electron located at one end. The individual buildup curves for the 40 nuclear spins are shown 
in F^. (III). The simulation parameters are listed in the following table (they are identical with the parameters used 
in [III) . To investigate the role of the nuclear dipolar coefficient dkj for the distribution of the polarization within 


Parameter 

Value 

Nuclei 


UJl 

36 MHz 

Temperature 

1 K 

Number of trajectories 

IcP 

CJl 

100 kHz 

|-^+|e —ni 

40 kHz 

|R+|e-n,-^l 

0 kHz 

A 

0 Hz 

(d) 

7.45 Hz ± 0.82 Hz 

Tie, T 2 e, tin, t 2 n 

10 ms, 10 ps, 10“ s, 1.25 ms 


the linear spin chain, we calculated the average polarization buildup and compared it to the buildup obtained when 
the diplar coefficient is set to half its value for all nuclei fc > 1 (fig. (|4|). In further simulations we either set all 
dipolar coefficients to half their value or ,in addition, we decreased the pseudo-secular interaction by a factor ^/2. 
These changes of the interaction parameters are in analogy to the set of simulations carried out by Hovav et al. in 
Ref. Ilf. In contrast to their observations, the biggest effect on the average polarization is already brought about 


by halving the dipolar interaction strength for all nuclei fc > 1, indicating the importance of this interaction for the 
diffusive transport. Any additional changes only introduced further small decreases of the average polarization. 

We then compared a simple diffusion model dp/dt = where p = p{x,t) is the polarization field and a 

constant source is assumed for the start of the linear spin chain. The average diffusion constant was obtained from 
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- normal 



Time fs^ Time 


FIG. 4. Polarization buildup in a linear chain of 40 coupled nuclear spins with an electron located at one end. The buildup 
curve for each individual nucleus is shown. On the right hand we show the average polarization for the spin chain for the cases 
that i) all dipolar constants are given by the spacing between the individual nuclei, ii) the dipolar coupling constants dkj are 
halved for bulk nuclei, iii) the dipolar coupling constants are halved for all nuclei and iv) the dipolar couplings are halved for 
all nuclei and in addition the pseudo-secular interaction strength between the electron and the first nuclei is reduced by \/2 


eq. ©. This diffusion model was solved for two different boundary conditions, either with an absorbing boundary 
condition or a reflective boundary condition. The time course of the polarization moving through the linear spin chain 
is compared in Fig. ([S]) with a contour plot derived from simulating the polarizationchanges within the linear spin 
chain using our method described in the main part of the paper. It is clear evident that the dynamics is reasonably 
well approximated by the simple diffusion model if a reflective boundary condition is used. The parameters for this 
simulation are provided in the following table. 


Parameter 

Value 

Nuclei 


Bo 

3.4 T 

Temperature 

1 K 

Number of trajectories 

10^ 

Wl 

50 kHz 

separation 

5 A ± 5% 

angle 9 of chain to Bo 

45° 

Tie, T 2 e, tin, t 2 n 

1 s, 10 /rs, 10“ s, 0.1 ms 



FIG. 5. Gomparison between a simple diffusion model and the results obtain by kMG simulations. The black contour lines 
visualise how nuclear spin polarization at various levels moves in time through a linear nuclear spin chain consisting of 30 
dipolar coupled spins during a SE DNP. These results are compared to a simple diffusion model that is based on a polarization 
source at the start of the spin chain and an average diffusion constant as described in eq. ©• The solid black line shows how 
polarization would move through the linear chain assuming a fully absorbing boundary at the end of the chain and the red line 
indicates the time course of the polarization if a reflective boundary condition is assumed. 






























